Isospin emission and flow at high baryon density: a test of the symmetry potential 



V. Giordanoi'2, M. Colonna\ M. Di Toroi'^.*^ V.Greco ^^^^ J. Rizzoi-^ 
^ Laboratori Nazionali del Sud, INFN, 
Via S. Sofia 62, 1-95123 Catama 
^ Physics and Astronomy Dept., 
University of Catania, Italy 
* email: ditoro@lns.infn.it 

High energy Heavy Ion Collisions (HIC) are studied in order to access nuclear matter properties at 
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implications. Using fully consistent transport simulations built on effective theories we test isospin 
observables ranging from nucleon/cluster emissions to collective flows (in particular the elliptic, 
squeeze out, part). The effects of the competition between stiffness and momentum dependence of 
the Symmetry Potential on the reaction dynamics are thoroughly analyzed. In this way we try to 
shed light on the controversial neutron/proton effective mass splitting at high baryon and isospin 
densities. New, more exclusive, experiments are suggested. 
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I. INTRODUCTION 

In the nuclear Equation of State (EoS) the symme- 
try energy Esym appears in the energy density e{p, pa) = 
e{p) + pEsymips/p)"^ + 0{p-i/ pY + •■, expressed in terms 
of total {p = pp + pn) and isospin (p^ = Pp — Pn) densi- 
ties. The symmetry term gets a kinetic contribution di- 
rectly from basic Pauli correlations and a potential part 
from the highly controversial isospin dependence of the 
effective interactions [H-ll]. Both at sub-saturation and 
supra-saturation densities, predictions based on the exist- 
ing many-body techniques diverge rather widely, see Q]. 
We remind that the knowledge of the EoS of asymmetric 
matter is very important at low densities (neutron skins, 
nuclear structure at the drip lines, neutron distillation 
in fragmentation, neutron star formation and crust..) as 
well as at high densities (transition to a deconfined phase, 
neutron star mass/radius, cooling, hybrid structure, for- 
mation of black holes...). We take advantage of recent 
opportunities in theory (development of rather reliable 
microscopic transport codes for HIC) and in experiments 
(availability of very asymmetric radioactive beams, im- 
proved possibility of measuring event- by- event correla- 
tions) to present results that could severely constrain the 
existing effective interaction models. We will discuss dis- 
sipative collisions at intermediate energies, i.e. in the 
range from IQQAMeV up to lAGeV. This will allow to 
probe the symmetry term above normal density. 

The reaction is simulated using transport codes 
based on effective mean field theories, with correla- 
tions/fluctuations included via hard nucleon-nucleon col- 
lisions, see [H, 0]. The dynamical effects of a different 
density dependence of the isovector part of the EoS are 
tested |7|]. Moreover a particular attention is devoted 
to the isospin effects on the momentum dependence of 
the symmetry potentials, i.e. to observables sensitive to 



a different neutron/proton effective mass in asymmetric 
matter. The problem of Momentum Dependence in the 
Isovector channel {Iso — MD) is still very controversial 
and it would be extremely important to get more definite 
experimental information, see the refs. [H, [1, [!]• Exotic 
beams at intermediate energies are of interest in order to 
have high momentum particles and to test regions of high 
baryon (isoscalar) and isospin (isovector) density during 
the reaction dynamics. 

The paper is arranged as follows. In Sect. 2 the self- 
consistent transport approach to the reaction dynamics 
is shortly described. Sect. 3 is devoted to a detailed anal- 
ysis of the used effective interactions, in particular of the 
isovector contributions to the local and non-local part 
of the mean field potentials. The latter, which corre- 
sponds to the Iso — MD case, leads to a splitting of the 
neutron-proton effective masses in asymmetric matter. 
Resuhs for Au -h^^^ Au are shown in Sect. 4, selecting 
the observables that would allow an independent study 
of the two isovector contributions. Conclusions and per- 
spectives are presented in Sect. 5. 



II. THE COLLISION DYNAMICS 

We perform ah initio collision simulations using the 
microscopic Stochastic Mean Field (SMF) model. It is 
based on mean field transport theory with correlations 
included via hard nucleon-nucleon (NN) collisions and 
with inclusion of stochastic forces acting on the mean 
phase-space trajectory tl, iS,, JTj . Stochasticity is essen- 
tial in order to allow the growth of dynamical instabilities 
leading to fragment production, as well as to obtain phys- 
ical widths of the observable distributions. Moreover it 
will allow to perform event-by-event correlation studies of 
great importance for the very complex reaction dynamics 
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in this energy range. 

The transport equation for the phase space distribu- 
tion function, with the Pauh blocking consistently evalu- 
ated, is integrated following a representation in terms of 
test particles of finite widths [l3| • A detailed descrip- 
tion of the procedure is given in ref. [2| • Our code (l2j 
has been extended by the introduction of momentum de- 
pendent mean fields (see next Section), which are rather 
important in this energy range. It has been also possible 
to improve the numerical accuracy while even reducing 
the computing times [Tsj . 

A parametrization of free nucleon-nucleon cross sec- 
tions is used, with isospin, energy and angular depen- 
dence [rij. Low energy NN collisions, mostly forbidden 
because of the Pauli blocking, have large cross sections 
and could induce spurious effects in the presence of some 
numerical inefficiency in the blocking procedure, due to 
the discretization of the phase space. In order to avoid 
such problems a cutoff value Ucut = 50to6 is used in our 
calculations. A parallel ensemble method is employed in 
the implementation of the collision term. 

For discussions of isospin dynamics in this energy 
regime it is essential to have a reliable procedure for frag- 
ment recognition, i.e. to identify free and clusterized nu- 
cleons. A coalescence procedure is applied in an event 
by event analysis at the "freeze-out" time, i.e. when the 
resulting fragments are well separated in space and in- 
teracting only via Coulomb forces. The fragments are 
formed using a phase-space proximity criterion. Two par- 
ticles are recognized to belong to the same cluster if they 
are sufficiently close in phase space. The procedure is ap- 
plied to several random samplings of N (total number of 
nucleons) test particles from the full ensemble distribu- 
tion at the freeze-out. The used coalescence parameters 
are dr — A.5fm, dp = 1.5/to~^, that reasonably repro- 
duce the charge distributions in the same energy range 
[isj . We have also seen that the results do not show ap- 
preciable changes for small variations of the coalescence 
parameters. Since the random sampling usually is not 
exactly preserving energy and momentum conservation, 
we apply an additional constraint to impose them. 



The kinetic term is: 

/d'^p 
-j^ [In {r, p)+fp (f, P}]^ (2) 

The terms £a and es of Eq.([T|) account for saturation 
properties, including the symmetry energy. 
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where p is the "local" density, pq the saturation value 
and /? = (iV - Z)/{N + Z) represents the "local" asym- 
metry parameter. The last term gives the momentum 
dependence and can be written as |17l |: 



8(C-h2z)^ 4(3C-4z), ^ 



bga 



(5) 



The momentum dependence is contained in the Irr' 
terms, which are integrals of the form 

Irr' = j dpdp'fr{r,p)fr'ir,p')g{p,p') , (6) 

with T — P, N , for protons and neutrons. Here /r(^jP) 
are the nucleon phase space distributions for protons and 
neutrons and the function g{p,p') = g[{p — p')'^] deter- 
mines the type of momentum dependence. A Skyrme-like 
momentum dependence is obtained when we use the sim- 
ple quadratic form g{p,p') = {p — p'Y' ■ A more general 
momentum dependence, in better agreement with phe- 
nomenological optical potentials, can be introduced by 
the function tl8n22i] 
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III. MOMENTUM DEPENDENCE OF THE 
FIELDS 



We adopt a generalized form of the effective interac- 
tions, which can be easily reduced to Skyrme-like forces, 
with momentum dependent terms also in the isovector 
channel ^, [H, [H, [l3|- The general structure of the 
isoscalar and isovector Momentum Dependent (MD) ef- 
fective fields is derived via an isospin asymmetric exten- 
sion of the Gale-Bertsch-DasGupta (GBD) force 18. - .22i] . 
which corresponds to a Yukawian non-locality. 

The energy density is parametrized as follows: 



e(r) = eH„(rO + sa (g(r)) + Eb (e(f)) + ec,z (rO) (1) 



For symmetric nuclear matter (p — 0) the energy den- 
sity Eqs. (|2I3I4I5I) reduces to the parametrization pro- 
posed by Gale, Bertsch and Das Gupta that for 
low densities and momen ta g ives the same results of the 
Skyrme-Gogny force, see |2l| . 

The EoS of symmetric matter is fixed by the pa- 
rameters A, B, C, a and A. We choose the values 
A = -111.3MeV, B = Ul.SMeV, C = -6A.5MeV, 
a — 7/6 and A — 1.5p^ (with Fermi momentum at 
normal density), which provide a rather Soft EoS with 
compressibility K ~ 215MeV and an isoscalar effective 
mass m* /m — 0.67 at saturation density po — 0.16/m~'^. 
We note that an indication for a Soft Eos (Isoscalar) even 
at higher densities comes also from heavy ion data (col- 
lective flows and meson production) at intermediate en- 
ergies [231] . 
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FIG. 1: Density dependence of the symmetry energies used in 
the simulations presented here: asy-soft (soUd) and asy-stiff 
(dashed). 



From the energy density one derives the mean field 
potentials as Ur{r,p) = Se/6fr- Thus the above energy 
density implies a momentum dependent mean field inter- 
action. The momentum dependence is isoscalar if, in the 
Eq.Jl]), the coefficients in front of Inp{r) and in front of 
the combination (/„„(r^ + equal (correspond- 

ing to the C = 8z case), but it can also get an isovector 
part, if they are different. 

The isovector momentum dependence implies different 
effective masses for protons and neutrons given as ^ — 
(1 + ?^^)~^ P = Pf.^> ^^'^ density. 
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FIG. 2: Upper Panels: Density dependence of the poten- 
tial symmetry energy (Solid Lines), in the Asy stiff choice. 
Dashed lines refer to local contributions, dot-dash lines to 
momentum-dependent ones, see text. Left: rrin > rrip 
parametrization; Right: < nip case. Lower Panels: cor- 
responding behavior of neutron/proton effective masses as a 
function of the density, for an asymmetry /? = 0.2. 

Here we probe the sensitivity of isospin observables 



to two essentially different density dependence of the 
symmetry energy around saturation: Asy-soft with a 
smooth behavior below saturation, and even a decrease 
at high densities, vs. Asy-stiff which always shows in- 
stead a rapid increase, roughly proportional to the den- 
sity [HiSIzI- In fact the two choices are well characterized 
by the slope parameter 

L = 3po{dEsym / dp) p=pg 

which is of the order of 10 — 20MeV in the Asysoft case, 
just from the kinetic contribution to the symmetry term, 
and around 70 — lOOMeV in the Asystiff parametriza- 
tion with also a repulsive potential part [H, 0]- In Fig-1 
we show the density dependence for these two typical 
choices. 

Table 1. Parameters xq, X3,z for opposite mass 
splittings, with the same Esym = i^MeV , for Asysoft 
(L=19MeV) and Asystiff (L=95MeV). 



Asysoft 








Mass splitting 






z{MeV) 


m; > to; 


1.111 


1.196 


-35.467 


TO* < to; 


5.212 


2.988 


44.971 


TO* =to; 


3.165 


2.094 


4.811 


Asystiff 








Mass splitting 


Xo 




z 


to; > to; 


-1.614 


-1.210 


-35.467 


to; < to; 


2.487 


0.582 


44.971 


TO„ = TOp 


0.440 


0.312 


4.811 



When we use momentum-dependent interactions we 
have also contributions to the symmetry energy from the 
non-local terms. The presence of an interplay between 
the C and z parameters in our form of the effective inter- 
action allows an independent study of the dynamical ef- 
fects of the stiffness of the symmetry term and of the neu- 
tron/proton effective mass splitting. This can be easily 
done just varying the xq, 2:3 and z parameters. In Table I 
we report the used set of parameters. The symmetry en- 
ergy at saturation is always fixed to Egymipo) = 33MeV. 

In Fig Owe plot the density dependence of the poten- 
tial part of the symmetry energy, in the Asystiff case, for 
the two choices of the n/p mass splitting (solid lines, up- 
per panels). We also separately report the contributions 
from the momentum-dependent, E{p,p), and the density 
dependent, E(p), part of the EoS, whose sum gives the 
total EPy^ . A change in the sign of the mass splitting is 
related to opposite behaviors of these two contributions, 
exactly like it happens in Skyrme-like forces, see sections 
(2.1-2.2) of ref.ilj. The lower panels show the density 
dependence of the corresponding mass splitting, for an 
asymmetry parameter /? = 0.2 (the ^^^Au asymmetry). 
In order to probe the mass splitting effects on the heavy 
ion dynamics we have chosen parametrizations that give 
almost opposite splittings at all densities. 

In Fig[3]we present the density dependence of the neu- 
tron/proton symmetry potentials, for the two stiffness of 
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FIG. 3: Density dependence of neutron(upper)-proton(lower) 
potentials for an asymmetry /3 = 0.2 for the Asysoft{left) 
and As j/sti// (right) choices. Dashed: No momentum depen- 
dence. Momentum dependent potentials at k — 2kF- solid 
lines for the > rrip case, dotted lines for the opposite 
ml < m* choice. 



the symmetry term, evaluated in the case without mo- 
mentum dependence (dashed lines) and in the momen- 
tum dependence {Iso — MD) case for the m* < m* (dot- 
ted) and the opposite m* > m* (soHd) choices. We see 
that the momentum dependence modifies the effect of the 
symmetry term stiffness on the nucleon potentials, with 
differences that become more appreciable with increas- 
ing nucleon momenta. From this figure we can already 
predict large effects of the effective mass splitting at high 
momenta. 

This is shown more explicitly in Fig|4] where we see 
the momentum dependence of the neutron-proton po- 
tentials at saturation density for the two mass split- 
ting choices, always for a "typical" /? = 0.2 asymmetry 
(I245'n^i97 Au...). The plot is for the Asysoft (left panel) 
and the Asystif f (right) symmetry term, and in fact it 
is not much different. Indeed we can see also from the 
previous Figl3]that at normal density the difference be- 
tween neutron and proton potentials is almost the same 
for the two asy-stiffness, even in the case of Iso — MD 
interactions.. 

The FigslS] and |4] suggest the presence of interesting 
Isoscalar and Isovector MD effects on the reaction dy- 
namics: 

• Isoscalar. In general the momentum dependence 
gives more attractive potentials at low momenta, 
p < Pf, and more repulsive at high p, p > pp. 
In the reaction dynamics we expect the more ener- 
getic nucleons to be fast emitted and to suffer less 
collisions [l8|. As a consequence we will have less 
stopping of the matter and less compression. The 
isoscalar EoS becomes stiffer when the momentum 
dependence is included. 



FIG. 4: Momentum dependence of neutron-proton potentials 
at saturation density and asymmetry P = 0.2, for the two 
sphtting choices ml^ < nip (dashed) and > rrip (solid). 
Left panel: Asysoft Iso — Eos. Right panel: Asystif f case 



Isovector. Isospin effects on the momentum depen- 
dence imply different slopes around pp for neutrons 
and protons, as clearly shown in FigUl and so the 
larger repulsion above is different. In the case 
m* < rn* the high momentum neutrons will see a 
more repulsive field with respect to the high-p pro- 
tons. The opposite will happen in the m* > m* 
case. The fast nucleon emission will be directly af- 
fected: in the to* < m* case we expect a larger n/p 
yield for nucleons emitted in central collisions and a 
larger neutron Squeeze — out (elliptic flow) in semi- 
central collisions in heavy ion reactions at interme- 
diate energies, in particular for high pt (transverse 
momentum) selections. In fact in the interacting, 
high density, early stage of the reaction dynamics 
the pressure is built from violent nucleon-nucleon 
collisions and the high pt particles will carry the 
maximal information on high density and momen- 
tum dependence of the symmetry potentials. The 
azimuthal distributions (elliptic flows) will be par- 
ticularly affected since particles mostly retain their 
high transverse momenta escaping along directions 
orthogonal to the reaction plane without suffering 
much rescattering processes. 

We will test those predictions also for n-rich vs. n- 
poor light ions, like i^H, ^He), easier to detect. 
Since, as already noted, the symmetry potentials 
are not very different in the Asystiff/ Asysoft choice 
for density range probed at intermediate energies 
(see the discussion in the next Section) , we can 
expect that the Mass — Splitting effect could be 
even larger than the one related to the different 
stiffness of the symmetry term. 

A final interesting point is about the crossing of 



FIG. 5: 197Au+197Au at 400AMeV (left) and 600AMeV 
(right), central collision. Longitudinal (solid lines) and 
Transversal (dashed lines) scaled rapidity distributions of 
Z=l-3 ions. 



the two predictions, with opposite mass-splitting 
choice, at high momenta, see FigH) Roughly this 
should happen around the Fermi momentum corre- 
sponding to the density reached in the compressed 
source of the fast emitted nucleons. This represents 
an independent way to check the maximum density 
reached during the collision. 



IV. RESULTS ON ^^'^ Au+^^'^ Au REACTIONS AT 
INTERMEDIATE ENERGIES 

Our Stochastic Mean Field (SMF) transport code has 
been implemented with I so — MD symmetry potentials, 
with a different {n,p) momentum dependence, as dis- 
cussed in detail in the previous section. This will allow to 
follow the dynamical effect of opposite n/p effective mass 
splitting while keepin g th e same density dependence of 
the symmetry energy [l7|. 
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We present here some results for ^^'^Ai 
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Au re- 
actions at AOO AM eV and 600AM eV [24[. For central 
collisions in the interacting zone we can reach baryon 
densities about 1.7 — 1.8po in a transient time of the or- 
der of 15-20 fm/c. The system is quickly expanding and 
the Freeze-Out time is around 50fm/c. At this time we 
have a dominant Coulomb interaction among the reac- 
tion products. All the results presented here refer to this 
time step. Secondary decays of excited primary frag- 
ments are not accounted for. In fact this will not affect 
too much the properties of nucleons and light ions at high 
transverse momenta mostly discussed in this work. 

For each beam energy and centrality we run 50 inde- 
pendent events. In order to show how global properties of 
the reaction are reproduced by the simulations, in FiglS] 
we report the longitudinal and transversal rapidity distri- 
butions of protons and light ions (Z=l-3) in central col- 
lisions at the two beam energies. We use the CM scaled 
rapidity yo = ycm/Uproj, vs. the projectile rapidity. A 
good parameter to evaluate the stopping power of the 
collision is the Vartl quantity, i.e. the ratio between the 
variances of the transverse and longitudinal rapidity dis- 
tributions, recently suggested by the FOPI collaboration, 
[l^Hl]. In our calculation we get values around 0.75, not 
much different at the two beam energies, comparable to 
data although no experimental trigger is included. We 





FIG. 6: 197Au+197Au at 400AMeV, central coUision.Isospin 
content of nucleon (left) and light ion (right) emissions vs. pt 
at midrapidity, | j/o |< 0.3, (upper) and kinetic energy (lower), 
for all rapidities, for the two nucleon mass splitting choices. 
Top Panels: Asysoft; Bottom Panels: Asystiff. 



note that this estimation of the stopping power appears 
not dependent on the isovector contributions to the EoS. 
This clearly shows the difficulty of the search for well 
measurable isospin effects in the reaction dynamics. Here 
we will focus our attention on the isospin content of fast 
emitted nucleons and light ions and on isospin collective 
flows. 



Isospin ratios of fast emitted particles 

In Figs l6l7l we present the (n/p) and ^H/^He yield ra- 
tios at freeze-out, for two choices of Asy-stiffness and 
Mass-splitting, vs. transverse momentum in a mid- 
rapidity selection (upper curves) and kinetic energy ( all 
rapidities, lower curves). In this way we can separate 
particle emissions from sources at different densities, as 
discussed in the previous Section. 

For both beam energies we clearly observe the opposite 
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FIG. 7: Same as in Fig.© for 197Au+197Au at 600AMeV, 
central collision. 



effect of the different mass splitting in the low and high 
momentum regions, as expected from FigH) E.g. in the 
771* < m* case the neutrons see a less repulsive potential 
at low momenta and a more repulsive one at high pt. The 
curves in the opposite mass-splitting show exactly the 
opposite behavior. We note some interesting features: 

i) The effect is almost not dependent on the stiffness of 
the symmetry term. At high pt, where particles mostly 
come from high density regions, the larger repulsion seen 
by neutrons in the Asy-stiff case, leading to an enhanced 
emission, is compensated by the larger Coulomb repul- 
sion in the remaining matter, favoring proton emission. 
On the other hand, at low pt, the sensitivity to the Asy- 
stiffness is lost due to the mixing of sources at different 
densities, also for central rapidities, during the radial ex- 
pansion. 

ii) The curves are crossing at pt — PprojectUe — 
2.13/m~^. The crossing nicely corresponds to the Fermi 
momentum of a source at baryon density p ~ l.Spoi 

iii) The results appear not very sensitive to the beam 
energy going from 400 to 600 AMeV, likely because the 
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FIG. 8: Proton (thick) and neutron (thin) V2 flows in a semi- 
central reaction Au+Au at 400AMeV. Top Panels: Rapidity 
dependence. Bottom Panels: Transverse momentum depen- 
dence at midrapidity, | j/o |< 0.3. Upper curves for m* > m*, 
lower curves for the opposite splitting < mj. Left: Asys- 
tiff. Right: Asysoft. 



reached maximum density is not much different in the 
two cases. The same is observed for the isospin flows 
analyzed in the following Subsection. 

We remark that all the effects discussed before should 
be also present for the ^H/^He yield ratios, more easily 
detected. Particularly interesting is the predicted large 
increase at high pt in the m* < m* choice. Some pre- 
liminary FOPf results seem to indicate this trend [2(|, 
but more data are needed. It is encouraging that we al- 
ready see a good Iso — MD dependence of rather inclu- 
sive nucleon/cluster emission data. In presence of a good 
statistics for the detection of high pt particles, a further 
selection at high azimuthal angles and central rapidities 
would certainly enhance the sensitivity to the momen- 
tum dependence of the Symmetry Potentials. This will 
introduce the discussion of isospin elliptic flows of the 
following Subsection. 
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FIG. 9: Same as for Fig[8] for the semi-central reaction 
Au+Au at eOOAMeV. 



Isospin flows 

Collective flows are very good candidates since they 
are expected to be very sensitive to the momentum de- 
pendence of the mean field, see [H [l^]- The flow ob- 
servables can be seen respectively as the first and second 
coefficients from the Fourier expansion of the azimuthal 
distribution: ^{y,Pt) ^No[l + Vicos{(t)) + 2V2Cos{2(t))], 

where pt = ypl+P^ is the transverse momentum and y 
the rapidity along beam direction. The transverse flow, 
Vi{y,pt) = provides information on the anisotropy 
of nucleon emission on the reaction plane. Very im- 
portant for the reaction dynamics is the elliptic flow, 

2 _ 2 

y2iy,Pt) = /" )■ The sign of V2 indicates the az- 
imuthal anisotropy of emission: on the reaction plane 
(V2 > 0) or out-of-plane {squeeze — out, V2 < 0) /Tj. 

Isospin effects on collective flows have been studied 
within the UrQMD transport model in order to probe 
the influence of the symmetry repulsion at high densities 
[28I [29} , here we focus the attention on the mass-splitting 
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FIG. 10: Transverse momentum dependence of the difference 
between proton and neutron V2 flows, at mid-rapidity, | |< 
0.3, in a semi-central reaction Au-|-Au at 400AMeV. Asystiff 
choice. 




FIG. 11: Proton (thick) and neutron (thin) V2 flows in a 
semi-central reaction Au-|-Au at 400AMeV, for equal (n.p)- 
effective masses. Top Panels: Rapidity dependence. Bottom 
Panels: Transverse momentum dependence at midrapidity. 
Left: Asystiff. Right: Asysoft. 



contributions. For the same Au-|-Au reactions, in a semi- 
central selection, we present in Figs. [51[H](Upper Panels) 
the rapidity dependence of (n/p) V2 for different choices 
of the Asy-stiffness and effective mass splitting. We ob- 
serve the relevance of the latter: at mid-rapidity the neu- 
tron squeeze-out is much larger in the m* < m* case in- 
dependently of the stiffness of the symmetry term. We 
note however that in the Asysoft case we see an inver- 
sion of the neutron/proton squeeze-out at mid-rapidity 
for the two effective mass-splittings. Good data seem to 
be suitable to disentangle Iso — MD potentials. 

The mass-splitting effect is large at high pt (Bottom 
Panels), again in a mid-rapidity selection, as expected for 
particle emitted from higher density regions. Here the 
results are also slightly depending on the Asy-stiffness, 
with large neutron squeeze — out effects in the Asystiff 
case. 
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FIG. 12: Rapidity dependence of Triton (thin) and ^He 
(thick) V2 fiows in a semi-central reaction Au+Au at 
400AMeV. Upper curves for m* > m*, lower curves for the 
opposite splitting < rrip. Left Panels: Asystiff. Right 
Panels: Asysoft. 




FIG. 13: Same as in Fig ll2l for a semi-central reaction Au-|-Au 
at eOOAMeV. 



A good sensitive observable seems to be the Difference 
of (p/n) elliptic flows, V2{p) — V2{n), shown in FiglTOlvs. 
transverse momentum at mid-rapidity, that is systemat- 
ically larger in the m* < m* case. 

In order to have a clear idea of the relevance of the (n,p) 
mass splitting on the fast nucleon emissions we present 
in FigHT] the neutron/proton elliptic flows for semicen- 
tral Au + Au collisions at AQQAMeV evaluated with the 
parametrizations giving to* = to* for the Au asym- 
metry /3 ~ 0.2 [so]- Now the isospin effects are only 
related to the different stiffness of the symmetry term at 
suprasaturation density. We see that, at variance with 
the mass-splitting results of FiglH the rapidity distribu- 
tions (top panels) are not much affected, with a slightly 
larger neutron squeeze-out in the Asystiff case. Consis- 
tently we see some difference in the transverse momen- 
tum dependence at mid-rapidity (bottom panels) only at 
very large pt- 

Due to the difficulties in measuring neutrons, we have 
analyzed the isospin sensitivity of light isobar flows, like 




FIG. 14: Transverse momentum dependence at mid-rapidity 
of '-^He (thick) and '"^ H (thin) V2 flows in a semi-central reac- 
tion Au-|-Au at 400AMeV. Upper curves for m* > m*, lower 
curves for the opposite splitting < nip. Left: Asystiff. 
Right: Asysoft. 



^Hvs. '^i/e and SO on. We still see effective mass splitting 
effects, although slightly reduced, as shown in Figs. [T^l 
[13] for the triton/^iJe case. 

As in the nucleon elliptic flow, at mid-rapidity the tri- 
ton squeeze-out is larger in the m* < to* case indepen- 
dently of the stiffness of the symmetry term. Again in 
the Asysoft case we see an inversion of the vs. ^He 
squeeze-out at mid-rapidity for the two choices of the 
mass-splitting. 

In Fig |14l we show also the the transverse momentum 
dependence of the ^He,^ H elliptic flows at mid-rapidity 
(I uo \< 0.3) at AOO AM eV. Some mass-splitting effects 
can be seen at high pt- This should be well observed 
in the flow difference. Unfortunately with the present 
number of events we cannot present a figure (like the 
FigCni of the {p, n) case) since we have too large error 
bars for the lack of statistics. 

We note that increasing the beam energy from 400 to 
600 AMeV the symmetry potential effects are not much 
changing, as already seen also for the isospin ratios of fast 
emitted particles discussed in the previous Subsection. In 
this respect an interesting positive result is coming from 
preliminary FOPI data on '^H vs. "^i/e flows, for Au-Au 
collisions at beam energies extended up to l.SAGeV. The 
triton V2 shows a larger squeeze-out at mid-rapidity (in 
a relatively high transverse momentum selection) |26|, 
clearly increasing with the beam energy, consistent with 
the reached higher baryon densities. 



V. CONCLUSION AND PERSPECTIVES 

The paper has been mainly devoted to the study 
of the effects of the momentum dependence of the in- 
medium effective nuclear interactions on heavy ion reac- 
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tions at intermediate energies. Apart the already known 
Isoscalar effects, like less stopping, less compression and 
an overall stiffer EoS, we have investigated in detail the 
Isovector contributions, that give raise to the poorly 
known neutron-proton effective mass splitting in asym- 
metric matter. We remind that this point is also of inter- 
est for relativistic bosonic models of the nuclear matter 
[H, 0, lil,, i3^] and the development of effective Lagrangi- 
nas for non-perturbative QCD related to the breaking of 
the chiral symmetry |33l - l35| . We have shown that the fast 
nucleon and light ions emissions in heavy ion collisions at 
intermediate energies are affected by the (n, p) mass split- 
ting. The isospin yield ratios {n/p,^ H/^He) and elliptic 
flows appear more sensitive to the Mass-splitting than 



to the stiffness of the Symmetry Term at high density, 
in particular for high transverse momentum selections. 
This opens the possibility of disentangling the two dis- 
tinct isospin effects on the reaction dynamics. New more 
accurate and exclusive data are needed. We like to men- 
tion the new measurements that will be p erformed at SIS- 
GSI by the ASYEOS Collaboration ^ and the new ex- 
periments planned at RIKEN- Tokyo and CSR-Lanzhou 
also with unstable, more neutron-rich, beams. 
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